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We study the growth of cosmological perturbations in the model of Lorentz-violating massive 
gravity. The Friedman equation in this model acquires an unconventional term due to the Lorentz- 
breaking condensates which has the equation of state w = —1/(37) with 7 being a free parameter 
taking values outside of the range [0,1/3]. Apart from the standard contributions, the perturba- 
tions above the Friedmann background contain an extra piece which is proportional to an arbitrary 
function i9(a; 1 ) of the space coordinates. This function appears as an integration constant and corre- 
sponds to a non-propagating scalar mode which may, however, become dynamical with the account 
of the higher-derivative corrections. For — 1 < 7 < and 7 = 1 the "anomalous" perturbations grow 
slower than the standard ones and thus the model is compatible with observations. Whether the 
model is experimentally acceptable at other values of 7 depends on the value of the function $(a;*) 
at the beginning of the radiation-dominated epoch. 
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I. INTRODUCTION 

The standard cosmological model is based on the as- 
sumption that the gravitational interaction is correctly 
described by the general relativity (GR) at scales com- 
parable to the horizon size. This model is quite success- 
ful in describing the bulk of the cosmological data: the 
anisotropies of the cosmic microwave background radia- 
tion, the primordial abundance of light elements and the 
structure formation in the early Universe. The quantita- 
tive agreement between the standard cosmological model 
and the observations has an ever-growing precision [1]. 
The question arises to which extent one should consider 
this agreement as a confirmation of the general relativity 
itself. 

In order to address this question an alternative model 
is needed whose predictions can be compared to those of 
GR. Such a model should coincide with GR at scales from 
~ 0.1 mm to the size of the solar system where GR has 
been tested directly. Therefore, the modifications of the 
gravitational interaction should occur in the infrared, at 
distances much larger than the size of the solar system. 

Interestingly, at these distances the predictions of GR 
actually do not agree with the observations; only after 
the introduction of the otherwise undetected dark matter 
and dark energy the agreement is achieved. The neces- 
sity for these new components is a major problem of the 
standard cosmology. Hence, in parallel with the direct 
searches for the dark components, the alternative mod- 
els of gravity should be explored which may eventually 
eliminate the need (or provide alternative candidates) for 
the dark matter and shed new light on the nature of the 
dark energy. 

It is a challenging problem to modify the gravitational 
interaction at large distances. An alternative model has 
to be theoretically consistent, i.e., free from ghosts and 



instabilities. In addition, it should be in agreement with 
the existing experimental data and should, ideally, pro- 
vide testable predictions for the future experiments. It is 
not obvious that such models exist, so it would be very 
important to construct an example. 

There have been several attempts made in this direc- 
tion [2, 3, 4, 5, 6, 7, 8, 9, 10]. In this paper we concen- 
trate on the massive gravity model [11] which is minimal 
in the sense that in does not contain new light prop- 
agating degrees of freedom as compared to the Einstein 
gravity. In this model the graviton acquires a mass due to 
the space-time dependent condensates of the four "Gold- 
stone" scalars (/>°(x), 4> 2 ( x )- The action of the model 
reads [12] 

S = J d 4 x^ [-M 2 Pl TZ + A 4 ^ (Z ij ) + £ mat tcr] , (1) 

where the first term is the standard Einstein-Hilbert ac- 
tion, £ ma tter stands for the minimally coupled ordinary 
matter and !F(Z l i) is a function of the derivatives of the 
four scalar fields (jP(x), (ftix) which depends on a single 
argument constructed as follows: 
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The constant 7 is a free parameter. The model possesses 
meaningful cosmological solutions for 7 > 1/3 and 7 < 
[13]. 

The flat space vacuum solution is [15] 



A 2 t, 
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The ground state is translationally invariant due to the 
derivative nature of the Goldstone coupling, but the 
Lorentz symmetry is spontaneously broken. The par- 
ticular dependence of the action (1) on the derivatives of 
the Goldstone fields through a single argument Z % i en- 
sures the non-pathological behavior of the perturbations 
about the vacuum solution [11], namely, the absence of 
ghosts and rapid instabilities. The tensor metric pertur- 
bations are, in general, massive with the mass determined 
by the first and the second derivatives of the function T . 
The low-energy spectrum consists of the two propagating 
tensor modes. The auxiliary scalars do not appear in the 
spectrum. 

An interesting and peculiar feature of this model is 
that the gravitational interaction between static sources 
is described by the standard Newton's law despite the 
non-zero mass of the graviton[16]. Due to this fea- 
ture the model passes the terrestrial and solar system 
tests even for graviton masses as large as (10 15 cm) -1 . 
Moreover, it admits the standard Fricdmann-Robertson- 
Walker cosmological solutions, the only trace of the Gold- 
stone scalars being two contributions to the energy den- 
sity which behave as a cosmological constant and as mat- 
ter with the equation of state w = — (37) -1 . The mas- 
sive gravitons may be created in cosmologically signifi- 
cant amount and may play a role of the dark matter [12]. 

Given that the model passes the most obvious con- 
straints one may wonder if it reproduces correctly more 
subtle parts of modern cosmology, in particular, the the- 
ory of structure formation. This is the question which we 
address in this paper. The answer is not obvious a priori 
since the vacuum in this model contains the condensates 
of the Goldstone fields whose perturbations mix with the 
matter density perturbations. 

We find that cosmological perturbations in the model 
(1) consist of two parts. The first part behaves identically 
to the perturbations in GR and is, therefore, compati- 
ble with observations. The second, "anomalous" part is 
proportional to an unknown function fi^x 1 ) of the space 
coordinates which arises as an integration constant. 

The appearance of this function reflects the existence 
of the non-dynamical mode with the dispersion relation 
lu' 2 = 0, much in common with the ghost condensate 
model [9]. In the model (1) the value of -d(x l ) is de- 
termined by the initial conditions. Note, however, that 
the action (1) is no more than the low-energy effective 
action. One should expect the higher-derivative correc- 
tions to eq. (1) to be present. In general, these corrections 
make a slowly changing dynamical variable. In this 
case the initial value of ^(x 1 ) may be determined by its 
evolution at the inflationary epoch. 

As we show below, the growth of the "anomalous" per- 
turbations depends on the value of 7. For — 1 < 7 < 
they grow slower than the standard ones, so that the 
latter dominate. At 7 = 1 the anomalous contributions 
to perturbations cancel out. Thus, at least in these two 
cases the perturbations behave in a standard way and 
the model (1) is consistent with the structure formation 



in the Universe. This is the main result of the paper. 

For other values of 7 (in particular, for the cosmo- 
logically interesting case 1/3 < 7 < 1 which corre- 
sponds to the equation of state of the condensate — 1 < 
w < —1/3) the anomalous perturbations grow at the 
radiation-dominated epoch, and may or may not grow 
faster than the standard perturbations at the matter- 
dominated stage, depending on the particular value of 
7 (see Sect. Ill for details). Whether this leads to a 
contradiction with observations depends on the function 
$(x 4 ) at the beginning of the radiation-dominated era. 
To address this question one needs to consider the higher- 
derivative corrections to the action (1). This will be done 
elsewhere. 

The paper is organized as follows. In Sect. II we review 
the cosmological solutions in massive gravity. In Sect. Ill 
we calculate the growth of perturbations in the FRW 
background. Sect. IV contains the summary of the results 
and their discussion. The details of the calculations are 
given in the Appendix. 

II. FRW SOLUTION 

In the absence of matter the model described by the 
action (1) admits the standard FRW solutions [13]. The 
spatially-flat homogeneous and isotropic ansatz [17] for 
the metric and Goldstone fields reads 

ds 2 = a 2 (77) (drj 2 — 8ijdx l dx^ , 
/ = A 2 ^), f = A¥, (3) 

which implies 

Z = ^/a 2 ^ 2 , 

where Z = —SijZ^/S and where prime denotes the 
derivative with respect to the conformal time rj. Assum- 
ing the ordinary matter is homogeneous and isotropic, 
the equations of motion (the Friedmann equation and 
the field equation for cf>°) take the form 

a 2 

K 2 = -^r(p m + P<f, + PA), (4) 

= d (a 3 " 1 ^ 1 - 1 / 2 ^) , (5) 

where TL = a' /a, so that Tl/a is the Hubble constant, 
and 2>Tz = S l ^dT jdZ l K The matter energy density p m 
has the standard form, while the two contributions of the 
Goldstone fields to the energy density read 

Pa = -A 4 T/2, P0 = -Z-i^ZTz. 

The first of these terms behaves like a cosmological con- 
stant, while the second contribution corresponds to mat- 
ter with the equation of state w = —1/(37). Thus, the 
Friedmann equation (4) reduces to the standard one, the 
only trace of the Goldstone field being the energy densi- 
ties p,), and pa- 
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For a given function T eq. (5) determines the depen- 
dence of the variable Z on the scale factor. For 7 > 1/3 
or 7 < this equation implies 

z i-i/2 7jCz ->0 as a -> 00. 

The case of interest (which occurs generically for an al- 
gebraic function T [13]) is when in this limit Z — > Zq = 
const, such that Ez(Zq) = 0. In this case the graviton 
mass remains finite at a — > 00 (see Sect. B for details). 

In a special case 7 = 1/3 the field equation (5) implies 
Z = const. Then Tz is not driven to zero by the cos- 
mological evolution and p^ behaves as the cosmological 
constant. The value of the total cosmological constant 
P4, + Pa is determined by the initial conditions. 

In what follows we neglect where possible the devia- 
tions from the point Z = Zq (note that in this point 
Ptj) = 0). We also assume that the value of pa is of the 
order of the present-day cosmological constant, and thus 
its contribution to the Friedmann equation at the epoch 
of structure formation is negligible. 

III. COSMOLOGICAL PERTURBATIONS 

As in the standard analysis of the cosmological pertur- 
bations [14], it is convenient to separate space and time 
components. The metric perturbations can be parame- 
terized in the following way, 

Sg 00 = 2a 2 ip, 

5g i0 = a 2 { Vl +d t B), (6) 
d gij = a 2 (2ip6 ij -d i F j -d j F i -2d i d j E + hi j ), 

where the vector perturbations Vi and Fi are transverse, 
while the tensor perturbation hij is transverse and trace- 
less. A similar parameterization can be used for the per- 
turbations of the Goldstone fields, 

5<f>° = A 2 e°, 

= A 2 (& + c^), (7) 

where £i is transverse. Finally, the perturbations of the 
ordinary matter are parameterized in the following way, 

<^u™ = (tiPm + Sp m ) UpU v - g^Spm - p m 5gp V 

+ (p m + Pm) {u u SUp + Up5u v ) , (8) 

where 5p m and 6p m are related by the matter equation 
of state and the perturbations of the velocity Sup are ex- 
pressed in terms of the scalar £ and the transverse vector 
£j as follows, 

8u t = a{d+ diC), 
Sua = Q><-P- 

Therefore, in total there is one tensor perturbation 
consisting of two components, four vectors Vi, F i: £j and 
Q, consisting of two components each, and 9 scalars ip, 



B, -0, E, £, C, 5p m and Sp m . One vector and two 
scalar perturbations are gauge degrees of freedom; they 
can be eliminated by imposing a gauge condition. As a 
consequence, there is only three gauge-invariant vector 
fields 

zui=Vi+F!, tn=ti-Fi, (9) 

and Q, and seven scalar gauge-invariant fields 

$ = ip-a- 1 [a(E' +B)]' , 

* = ip + H(E' + B), 

E = £ — E, 

B = B + E'-f/cj)', 

S P = [8 Pm - P ' m {E' + B))/p m , 

S P = [S Pm -p' m (E' + B)]/p m , 

S< = (-(E' + B). (10) 

The tensor perturbation hij is also gauge invariant. 

A. The tensor perturbations 

The equation for the tensor perturbation is 

h'lj + ZHh^ - dfh i:j + afmlhij = 0, (11) 

where mi is the graviton mass which has the scale 
TO2 oc A 2 /Mpi and whose precise expression in terms of 
the function T(Z % i} and its derivatives is given in the 
Appendix B. 

This equation is identical to the equation for a free 
massive scalar field in the FRW background. If the mass 
of the graviton is larger than the Hubble constant, 3> 
H/a, which we assume to be the case in what follows, 
eq. (11) describes massive gravitational waves with the 
amplitude scaling like hij « a" 1 and hij oc a~ 3 / 2 in the 
relativistic and non-relativistic limits, respectively 

B. The vector perturbations 

In the longitudinal gauge the three equations describ- 
ing vector perturbations are 

d)wi - 2a 2 p m M- 2 {\ + w)( t = 0, (12) 
137^ + 2TLwi — a 2 m\<Ji = 0, (13) 
m\d)<n = 0, (14) 

where w is the parameter entering the equation of state 
of the ordinary matter, p m — wp m . 

The first of these equations allows to express £ l in 
terms of Wi, while the third equation gives cr.; = 0. There- 
fore, the only non-trivial equation is eq. (13). It differs 
from the conventional one by the term proportional to 
the graviton mass ra 2 . which cancels at Cj = 0. Thus, 
this equation is the conventional one and describes a field 
with the amplitude decreasing as Wi oc a~ 2 . 
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C. The scalar perturbations 

The perturbation S p can be expressed in terms of 5 P 
by means of the matter equation of state. In the case of 
the adiabatic perturbations one has 



S p = —S„, 
w 

where c s is the sound velocity (c 2 = w for the ideal fluid). 
The behavior of the remaining six perturbations is gov- 
erned by the following equations, 



Here m, are the graviton mass parameters defined in the 
Appendix B. We assume that they are of the order of 
A 2 /M P1 . 



This system of equations can be solved as follows. At 
mf 7^ eq. (20) can be used to express <f>" /4>'B + B' + $. 
Then eq. (19) becomes a closed equation for 3 
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o = -2a;* 

2 Pm 



m {H$ + *') + a 2 m 2 (d]E 



-a' 



1v1 pi 



a m 



h-B + B' + $ 



= -2*" - 2$ (H 2 + 2H') + d 2 (* - $) 



-2^(2* + $) 



a 



A/2 



m 2 d 2 ~ 



+m\ $ + j-B + B' - */7 



3* 



(15) 
(16) 

) 

(17) 



(18) 



The solution of this equation which docs not grow at 
spatial infinity reads 



(21) 



where $ (x 1 } is an integration constant depending only 
on spatial coordinates. Then eqs. (16), (17) and (20) can 
be used to express S^, S p and B in terms of $ and * while 
eq. (15) reads 



o = d Q 



,4^,2 



37 ( $ + B' 



-B 



-(3* + a?s))], 

ml I */7- ^ttB-B' -* 
V 

+ (ml - m|) 9|3. 



(19) 



(20) 



(22) 



With the account of all these relations, the remaining 
eq. (18) becomes a closed inhomogeneous equation for 



a 



w\ 1 

77T W + - 



H 2 



H 2 
1 

1 



2^ 



(1 + 



3c: 



2^ 

H 2 



H 1 



(23) 



Once the solution to this equation is found, the other 
variables are determined by eqs. (22), (16), (17) and (20). 
In particular, one finds 



S P = 



2M 2 t 



2M 2 
a 2 p m 



m z 1 



m 2 ) a 177 " 3 ?? 
d 



(24) 



where * is a solution to eq. (23). 

The conventional cosmological perturbations are re- 
covered by setting the graviton masses to zero, mf = 0. 



In this case eq. (15) gives $ — \t = which implies 
i?(x 4 ) = (cf. eqs. (21) and (22)). Then both cq. (23) and 
eq. (24) reduce to the standard equations describing cos- 
mological perturbations in the Einstein theory. Note that 
the value of ■d is determined essentially by the initial con- 
ditions. Setting $ = would eliminate the i?-dcpcndent 
terms in eqs. (23) and (24) and bring these equations to 
the conventional form even in the case m\ ^ 0. 



In the case of matter perturbations in a matter- 
dominated Universe eq. (23) reduces to the following 
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equation, 

<9 2 * 7 9* fl 



+ J_?£. + ±_1 o V7-30 = Oj 



9a 2 2a da \ 7 



which differs from the standard case by the presence of 
the inhomogeneous term proportional to •&. The solution 
to this equation reads 



2 + 3 7 



where Ci(x z ) are the integration constants. Substituting 
this solution into eq. (24) one finds the density contrast 



S P = 



Po 



^^-■(^.^J. (25) 



where po is the energy density of matter at present. The 
first two terms in this equation are precisely the ones 
which appear in the standard Einstein theory, the sec- 
ond term describing the linear growth of the perturba- 
tions, S p oc a. The difference with the conventional case 
consists in the third term on the r.h.s of cq. (25). The 
perturbations corresponding to this term grow propor- 
tionally to a 1 / 7 . For 7 > 1 or 7 < these "anomalous" 
perturbations grow slower than the standard ones. 

In the radiation epoch the situation is similar. For 
a relativistic fluid one has w m = 1/3, so that eq. (23) 
reduces to the following one, 



a 2 * 




4d* 
a da 

- 1 



m 2 p A 

Pr 
Pr 



i 1 / 7 " 3 ^, (26) 



where p r is the energy density of radiation at present. 
For the generic value of 7 the solution to this equation 
is cumbersome. For simplicity let us concentrate on the 
modes which are much smaller than the Hubble scale, 
k 2 ^> Ti. 2 . The density contrast calculated according to 
eq. (24) has the standard oscillating piece and the extra 
part proportional to 



S p ~ c\(x l ) siny + C2{x l ) cosy + 27 



k 2 M 2 pl 



-V 1+lh + / dxx 1+1 ^sin{y-x) 
Jo 



(l/7-l)/2 



0, (27) 



where y = r\kj\f^> is proportional to the scale factor, 
while Ci^x 1 ) are two integration constants. As one may 
see from this expression, for — 1 < 7 < the ^-dependent 
contribution to the density contrast decays with the scale 
factor so that only the standard contribution remains. 



Thus, in this range of 7 the perturbations behave just 
as predicted by general relativity in both matter and 
radiation-dominated epochs. 

Another case of interest is 7 = 1. This case is spe- 
cial because at 7 = 1 the a-dependence of the last term 
in eq. (26) disappears. In fact, one may show that in 
this case the dependence on $ cancels out in the density 
contrast, so that only the standard part of perturbations 
remains. 

At other values of 7 the independent contributions to 
perturbations grow in the radiation-dominated Universe. 



IV. DISCUSSION 

To summarize, in the model of massive gravity de- 
scribed by the action (1) the cosmological perturbations 
contain two parts, the "normal" and the "anomalous" 
one. The first, normal part has the behavior identical 
to that found in the conventional general relativity. It 
is therefore in agreement with observations to the same 
extent as the latter. In particular, the "normal" part of 
the perturbations can describe successfully at least the 
linear stage of the structure formation. In fact, the pre- 
dictions of the general relativity and the model (1) with 
only the "normal" perturbations present completely co- 
incide in the linear regime, so that the two models are 
indistinguishable in this respect. 

The second, "anomalous" part of perturbations is spe- 
cific to the model of massive gravity with the action (1). 
These perturbations originate from the condensates of 
the scalar fields present in the model. This contribu- 
tion to perturbations depends linearly on the unknown 
function of space coordinates i)(x l ) which enters as an in- 
tegration constant. The value of this function cannot be 
determined within the model (1) and has to be specified 
as an initial conditions. 

The behavior of the "anomalous" contributions to per- 
turbations at different stages of the evolution of the 
Universe depends on the value of 7. At the matter- 
dominated stage the anomalous perturbations grow not 
faster than the standard ones for 7 > 1 and 7 < 0. In 
the radiation-dominated epoch this occurs at — 1 < 7 < 
and 7=1. Thus, at —1 < 7 < and 7 = 1 the normal 
perturbations dominate at both radiation and matter- 
dominated stages. 

The appearance of the time-independent arbitrary 
function is not surprising. The same function ^(x 1 ) also 
enters the expression for the gravitational potential of 
an isolated massive body [12]. Its origin may be traced 
back to the existence of the scalar mode with the disper- 
sion relation lo 2 = [11]. This mode is not dynamical in 
the model (1). However, the action (1) is the low-energy 
effective action, so one should expect the corrections con- 
taining higher-derivative terms to be present. In general, 
these corrections make 1? a dynamical variable with the 
dispersion relation ui 2 = ap 4 , where a is a small coef- 
ficient. Therefore, $ becomes a slowly varying function 



G 



of time. The slow evolution may drive $ to a particu- 
lar value at the inflationary epoch and thus prepare the 
initial conditions for the radiation-dominated stage. If 
this initial value of i? is small, then the growth of the 
"anomalous" part of perturbations may become irrele- 
vant and corresponding values of 7 phcnomenologically 
acceptable. This question will be considered elsewhere. 
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For the cosmological solutions given by cq. (3), X = 
a -V 2 , V i = and W lj = -a- 2 5 lj . The 14 field equa- 
tions reduce to the following three relations (prime stands 
for the derivative with respect to 77), 

a 2 

= V72" (Pm + P<t> + Pa) , 
lvl P i 

2H' + H 2 = (Pm+P*+ Pa), 

lvl P i 

= d (a?F x xV*) , 

which can be solved for any given function T (a, X). The 
following notations were used 

9<s> = A±XT X , P<p = A 4 WT W , 
Pa = -A 4 T/2, p\ = A 4 T/2. 



APPENDIX 

In these Appendix we provide some intermediate for- 
mulae skipped in the main text. In Appendix A we cal- 
culate the field equations for the standard FRW solution 
(3). In Appendix B we introduce the mass parameters 
and the relations between them. In Appendix we calcu- 
late the gauge- invarant linearized field equations. 



For the model (1) characterized by a function T depend- 
ing only on Z %3 one has — — (37) p^. 



APPENDIX B: THE MASS PARAMETERS 

There are five mass parameters m 2 , i = . . . 4 defined 
by the following relations 



APPENDIX A: THE BACKGROUND 

For X, V 1 and W %3 given by eq. (2), the following 
action describes a massive gravitational field 

S = J d 4 xv^g [-M 2 P {R + A 4 T (X, V\W ij , ...)], 

characterized by the ten Einstein equations Q^ v = 
M~ 2 (T™ + T* v ) and four Goldstone equations 



-99 



a/3 



8X 



V 3 V % V 3 
-—d a {<f>\ + + 5l—d a( j? 



Here T™ is the energy-momentum tensor of the ordinary 
matter which has the standard form (8) while the energy- 
momentum tensor of the Goldstone field is given by 

7*, = A 4 {-\g»,? (^<W° + d^W) 



+ 



V 3 

drfduP - — {d^d^° + dp<jPd v <t> 1 ) 



dT dT o 
1 —du(f> du4> 



V l V 3 

+^ r d^°d^° 



8W i3 dX 



A 4 

1 2 - — — 

~~ M 2 
1u pi 

2A 4 



M 2 
ivi pl 

2A 4 



A 4 



3 ~ M 2 PI 
A 4 



(XT X + 2X 2 T X x) , 



-XT X - WT W + ^XWT vv ^j , 
(WT W - 2W 2 T WW2 ) , 
(WT W + 2W 2 T WW1 ), 



Af 2 
IU Pi 



[XT X + 2XWTxw) , 



where W = —l/SSijW" 13 and where the first and sec- 
ond derivatives of the function T (X, V % , W lJ , . . .) are 
denoted as follows 



dT 



— J~x, 



OX 
d 2 T 

dV l V 3 
dT 

dW 13 
d 2 T 
dXW 13 
d 2 T 
dW i3 W kl 



dX 

d 2 T 
— 2 = -rxx, 



= Tvv^ij, 

= Tw^ij, 

= TxwSij, 

= TwWlSijSkl + TwW2 (SikSjl + SuSjk) 
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It is possible to relate the masses to 2 , and m\ by using the 
Goldstone equation for the background. This last gives 



ml ( - H ) = mm\. 



For the model characterized by the function T 
T (Z 1 ^, it is straightforward to show that 



A 4 
Ml 



7 [3(1 - 2 1 )ZT Z 



+6 1 Z 2 (3F ZZ1 + 2T ZZ2 )1 
2A 4 

_ r (3 7 -l)^, 
lvl P i 

2A 4 , 2 
—3- (ZJ^z - 2Z 2 ^j Z2 ) , 



A 4 
A 4 



[ZTz + 2Z 2 Tzz\) , 
7 [Z.F Z + 2Z 2 (3^zzi + 2^zz 2 )] 



where 



QT 
dZv 
2 T 

dZ l 3Z kl 



^FzziSijSki + Tzzi (SikS-ji + SaSjk) 



For this particular case, — 3jp^, 
relations hold 



and the following 



stein tensor is given by 

sg oa = 2%*- mi/, 

6Q i0 = 0, [2 (Hip + lb') + (2W + H 2 ) B] + ^0 2 w t 

+ Vl {2H 1 + H 2 ) , 
SQy = -5 i:j [2 (<p + i/>) (2H 2 + W) - 2H (2ip + <p)' 

-2ii)" + o 2 (*-*)- 6 & + ip) (n 2 + h')] 

+didj [vP - $ - 2 (pa + U 2 ) E] - Htiy 
+ \ {dfhij - h'Q + (2H' + H 2 ) ha 

+H (djZUj + d-jZUi) + - (dizu'j + 3 vo'/) 
-(H 2 + 2H') (diFj+djFi). 

The linearized energy-momentum tensor of the ordinary 
matter fields is given by eq. (8), while the linearized en- 
ergy momentum tensor of the Goldstone fields takes the 
form 



5% 



nil 



a 2 M 2 Pl [(2 P4> + 2p A -m 2 )ip + n$£ u 
+ (/50+P0 + m 2 ) (d 2 Z + 3^)] , 



81% = 



M 2 Pl 



'4> + PA 

M 2 



{Vi + d,B) 



+-^(vi + diB + & + a&) 



a 2 M 2 



pi i 



a 2 M 2 Pl 



-m 2 (drfi + <U\ + 2<),i),0 



m 2 <% (3^ + 2 Z) + 5 l3 m\ {tp - £°70') 

\ ml - PJ ^if) own - w - W 

20,0,1-: • //,,:■ . 



ml = 37 (m\ - m\/2) , 
ml = 7 (3m 2 - rnj) , 
m\ = 2(3 7 -l)p . 



These relations allow to write the linearized Einstein 
field equations for the massive gravitational field. They 
consist of one tensor equation 



= h'fj - dfh i:j + 2Wiy + a 2 m 2 %, 



two vector equations 



APPENDIX C: THE LINEARIZED EQUATIONS 



= + 2TLvji — a 2 m\(Ji 



Small perturbations on the cosmological background 
can be expressed trough eq. (6-7). The tensor field 
and the vector field Q are gauge-invariant, while the 
other two gauge-invariant vector fields and seven gauge- 
invariant scalar fields are introduced trough eq. (9-10). 
With all these definitions, the linearized part of the Ein- 



a 2 2 Wi 



-2 



M 2 Pl 



(1 + UJ m ) Cii 



8 



and four scalar equations 
= -2d 2 * + m (H® + *') - a 2 ml 

(d 2 E + 3* 



M 



Pl 



2 P 



M 2 Pl 



dp, 



= 2d, (*' + H$) + a 2 d t 
P4> + P4, „, p 



M 2 



M 2 



(l+w U! )S ( 



= didj ($ - * + a 2 mjZ) , 



-2#" - 2$ (H 2 + 2U 1 ) + d 2 (* - $) 



-2H (2* + $)' + a 2 P 



M 2 p 
iu pi 



2 2 ^2 — 

Or, — a m,o, a 



+a 2 m 2 ^$ + ^B + B'J + a 2 (to 2 ~ 3m 2 ) 
The linearized Goldstonc equations give one vector equa- 



4 U? + 2^ 



tion 

= a~ 4 a 

and two scalar equations 
= do 



mld 2 (Ji, 



to 2 (^-B + B' + ^ ^m 2 (3^ + d 2 E 



d t d { a 4 

- a 4 a 



M 



pi 



ml - to 2 ) d 2 ^ 



AF 2 

1V1 PI 

2 , P<t>+P<l> 

4 M 2 



d> 



[1] D. N. Spergel et al. (WMAP), Astrophys. J. Suppl. 170, 

377 (2007), astro-ph/0603449. 
[2] M. Milgrom, Astrophys. J. 270, 371 (1983). 
[3] J. D. Bekenstein, Phys. Rev. D70, 083509 (2004), astro- 

ph/0403694. 

[4] R. Gregory, V. A. Rubakov, and S. M. Sibiryakov, Phys. 

Rev. Lett. 84, 5928 (2000), hep-th/0002072. 
[5] G. R. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. 

B485, 208 (2000), hep-th/0005016. 
[6] S. M. Carroll, V. Duvvuri, M. Trodden, and M. S. Turner, 

Phys. Rev. D70, 043528 (2004), astro-ph/0306438. 
[7] I. I. Kogan, S. Mouslopoulos, and A. Papazoglou, Phys. 

Lett. B501, 140 (2001), hcp-th/0011141. 
[8] T. Damour and I. I. Kogan, Phys. Rev. D66, 104024 

(2002), hep-th/0206042. 
[9] N. Arkani-Hamed, H.-C. Cheng, M. A. Luty, and 

S. Mukohyama, JHEP 05, 074 (2004), hep-th/0312099. 
[10] V. A. Rubakov (2004), hcp-th/0407104. 



[11] S. L. Dubovsky, JHEP 10, 076 (2004), hep-th/0409124. 

[12] S. L. Dubovsky, P. G. Tinyakov, and I. I. Tkachev, Phys. 
Rev. Lett. 94, 181102 (2005), hep-th/0411158. 

[13] S. L. Dubovsky, P. G. Tinyakov, and I. I. Tkachev, Phys. 
Rev. D72, 084011 (2005), hep-th/0504067. 

[14] V. F. Mukhanov, H. A. Feldman, and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992). 

[15] To be precise, the vacuum solution reads 4>° = at, 4>* = 
bx 1 , where a and b are two constants whose value is set by 
the requirement that the energy-momentum of the scalar 
fields in the Minkowski space vanishes. We assume in 
what follows that the function T is such that a — b = A 2 . 

[16] This feature is obviously a consequence of the breakdown 
of the Lorentz invariance. 

[17] This ansatz cannot be directly generalized to the cases of 
open and closed Universe [11]. It is not known whether 
the model admits cosmological solutions of these types. 



